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ABSTRACT : 

The  results  of  an  analog  computer  study  of  hypersonic  laminar  heat 
transfer  including  the  effects  of  pressure  gradient,  mass  transfer, 
velocity  slip,  Prandtl  number,  and  wall  temperature  are  presented. 
Reynolds  analogy  factors  and  adiabatic  wall  enthalpy  results  are  included 
for  all  combinations  of  the  aforementioned  effects.   In  addition,  an 
auxiliary  correlation  is  presented  which  simplifies  the  specification  of 
the  Reynolds  analogy  factor  with  velocity  slip.  All  of  the  computations 
presented  in  this  paper  were  performed  on  the  EAI  580  analog  computer  in 
the  Department  of  Aeronautics . 
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NOMENCLATURE 
(3  =  Pressure  gradient  parameter  in  equation  (l) 

C  =  Skin  friction  coefficient     p,  r—  \l-x   p  u 

C  =  Specific  heat  at  constant  pressure 

r^peU   y 
71  =  Independent  variable     f  -=-  dy 

Jk    Jop- 


e 


f *  =  Normalized  velocity     U/U 


e 


g  =  Normalized  enthalpy     H/H 

H  =  Stagnation  enthalpy     C  T  +  IT/2 

j  =  Dimensionality  parameter  0  for  two  dimensions,  1  for  axisymmetric 

k  =  Coefficient  of  thermal  conductivity 

(j,  =  Coefficient  of  viscosity 

Pr  =  Prandtl  number     llC  /k 

P 

o 
q  =  Heat  transfer  rate  in  BTU/ft  -sec 

r  =  Local  body  radius 

p   =  Fluid  density 

S  =  Reynolds  analogy  fact  defined  in  equations  (6a)  and  (7a) 

S*  =  Correlation  defined  in  equation  (8) 

St  =  Stanton  number  defined  in  equation  (k) 

U  =  Streamwise  velocity  component 

x  =  Streamwise  coordinate 

y  =  Lateral  coordinate 

£  =  Independent  variable     J   p  |i  U  r  °dx 


IV 


SUBSCRIPTS 
aw  =  Adiabatic  wall  conditions 
e  =  Local  external  conditions 
w  =  Wall  conditions 

SUPERSCRIPTS 


i   _ 


Differentiation  with  respect  to  T) 


v 


INTRODUCTION 
Under  the   assumption  of  self-similarity,  the  equations  of  momentum 
and  energy  for  a  laminar  "boundary  layer  in  a  locally  hypersonic  flow  can 
be  written  as 

f"»    +  ff"   +  p(g  -  f'2)   =  o  (1) 

g"    +  Pr  fg'    +  2(Pr   -  1)    (ff")'    =  0  (2) 

subject  to  the  boundary  conditions 

f     =  f 

w 

@  T]  =  0  f'    =  f  * 

1  w 


g     =gworg'=  g'  (3) 


@     7]     =     » 


f '  -  1 
ff   -  1 


Equations  (l)  and  (2)  describe  the  flow  of  a  constant  property,  perfect  gas 
whose  viscosity  varies  linearly  with  temperature.   In  addition,  the 
restriction  that  the  pressure  gradient  parameter,  p,  be  a  constant,  implies 
that  the  pressure  varies  linearly  with  streamwise  distance. 

The  results  of  an  analog  computer  study  of  hypersonic  laminar  heat 
transfer  involving  some  fifteen  hundred  solutions  of  equations  (l)  and  (2) 
are  presented.  These  calculations  provide  necessary  data  for  making 
estimates  of  the  convective  heat  transfer  rates  utilizing  Reynolds  analogy 
in  situations  involving  mass  transfer,  velocity  slip  and  strong  pressure 
gradients . 

A  comparison  of  the  analog  solution  of  the  Falkner-Skan  equations  with 
the  numerical  values  of  Cohen  and  Reshotko  [l]  indicated  an  error  of  less 


than  one  percent.  Because  of  the  computational  complexities  introduced 
by  the  additional  nonlinear it ies  of  equations  (l)  and  (2),  the  data 
presented  is  estimated  to  be  accurate  to  within  two  percent. 

DISCUSSION 
Based  on  Reynolds  analogy,  the  following  equations  relate  the  local 
convective  heat  transfer  rate  to  the  wall  shear  and  the  driving  enthalpy 
difference 

%   =  PeUeHe  <«*  "  *J   St  <*) 

St  =  Cf/2S  (5) 

In  this  formulation,  g   represents  the  recovery  factor  which  was  calcu- 
lated for  all  combinations  of  Pr,  f  ,  and  (3  in  the  range  of  values 
studied.  These  results,  presented  in  Fig.  1,  reinforce  the  classical 
result  that  the  laminar  recovery  factor  varies  as  the  square  root  of  the 
Prandtl  number  for  zero  mass  transfer.  However  with  blowing  or  suction 
no  such  simple  Prandtl  number  correlation  is  possible. 

When  calculating  heat  transfer  using  Reynolds  analogy,  the  Reynolds 
analogy  factor,  S,  which  exhibits  sensitivity  to  all  the  parameters  of 
the  problem,  must  also  be  specified.  Therefore,  calculated  values  of 
Reynolds  analogy  factors  are  presented  for  all  combinations  of  g,  =  1.1, 

1.0,  .5,  0;  Pr  =  1.25,  1.0,  .72,  .5;  f  =  0,  .25,  -.25;  P  =  -.2  -  1.0 

w 

for  f '   =0.   In  these  cases  we  can  compute  the  Reynolds  analogy  factor 
from  the  calculated  data  with  the  equation 


B-Prf  (gw  ,"  gg)  (6a) 


w 

w  w 


where,  in  conjunction  with  equation  (h) ,  we  have  utilized  Fourier's  law 
at  the  wall 

dT 


%  -  "kw  §7  i  <7a) 


W 

In  addition  to  these  situations,  a  case  of  large  velocity  slip  (f  =  .1) 

W 

is  included  for  Pr  =   .72  with  all  other  parameters  varied  as  before.      In 
this  case,  equation   (6a)   must  be  modified  to  account  for  velocity  slip  and 
takes  the  form 


(g       -  e  )   f"     Pr 
e  _  aw       *w'         w ,(■    n 

S  "  g»      +  2f   f"      (Pr  -  1)  (6b) 


W  WW 


Which  reduces  to  equation  (6a)  when  f '  =  0  or  for  a  Prandtl  number  of  one 

with  an  arbitrary  f  .  To  obtain  equation  (6b)  we  have  used  the  definition 

of  H  at  the  wall  and  Maslen's  [2]  results  that  the  heat  transfer  at  the 
wall  must  now  be  written 

%  -  -  K 1 1  +^A  1 1  )  '                  <7b) 

w         w 


The  results  of  the  computations  of  the  Reynolds  analogy  factor  without 
velocity  slip  are  plotted  in  Figs.  2a  -  21.  The  data  deviates  from  classi- 
cal results  for  nonhypersonic  boundary  layers  and  the  range  of  values  which 
the  Reynolds  analogy  factor  encompasses  illustrates  the  potential  uncer- 
tainty in  the  calculation  of  heat  transfer  without  knowledge  of  the 
dependence  of  S  on  the  particular  parameters  of  the  problem. 

The  data  was  processed  to  evaluate  various  Pr  -  (3  correlations  similar 
to  that  of  Tifford  and  Chu  [3]  in  an  effort  to  coalesce  the  dependence  of 
S  on  the  variables  of  the  problem.  No  satisfactory  general  correlation 
could  be  found  and  specific  information  on  the  appropriate  Reynolds  analogy 
factor  must  be  obtained. 


For  the  cases  involving  velocity  slip  without  mass  transfer  and  with 
Prandtl  number  different  than  unity,  the  data  was  processed  with  an 
auxiliary  function 

g   -  g, 

S*  =  f"   Pr  (  aw  ,    +  2f'  )  (8) 

w    v   g'        w'  v  ' 

&  w 

The  results  of  these  computations  are  presented  in  Fig.  3-  When  these 

plots  of  S*  are  compared  to  the  corresponding  plots  of  Reynolds  analogy 

factor  without  velocity  slip  for  the  same  Prandtl  number  and  wall  enthalpy, 

it  is  seen  that  the  curves  agree  to  within  a  few  percent.  Therefore,  from 

knowledge  of  S*,  the  Reynolds  analogy  factor  with  velocity  slip  may  be 

obtained  from  equations  (6)  and  (8) 

2f ' 
1  =       1         +     w     Pr  -  1  r.x 

S   S*  -  2f*   f1"1   Pr   g   -  g   *   Pr  ^ 

w    w     °aw    w 

where  values  for  ?*   are  taken  from  corresponding  values  of  S  for  the 
appropriate  Prandtl  number  and  wall  enthalpy  but  without  velocity  slip. 
Note  that  for  a  Prandtl  number  of  one,  the  correlation  is  both  inaccurate 
and  unnecessary  since  equation  (6b)  says  that  S  is  unaffected  by  velocity 
slip  in  this  case. 

CONCLUSIONS 
The  wide  range  of  values  which  the  Reynolds  analogy  factor  assumes  in 
hypersonic  flow  with  mass  transfer  and  strong  pressure  gradients  forces 
one  to  obtain  data  on  the  specific  situation  of  interest.  The  analog 
computer  has  yielded  fast,  accurate  solutions  of  the  equations  to  provide 
this  data  and  investigate  all  important  effects.   In  addition  an  approxi- 
mate correlation  has  been  found  which  reduces  the  need  for  detailed  inform- 
ation on  Reynolds  analogy  factors  in  a  flow  with  velocity  slip. 
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Fig.  2a  Reynolds  Analogy  Factor  vs.  0  for  g^  =  0 
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Fig.  2  c  Reynolds  Analogy  Factor  vs.  p  for  g  =  0 
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Fig.  2d  Reynolds  Analogy  Factor  vs.  p  for  g  =  .5 
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Fig.  2e  Reynolds  Analogy  Factor  vs.  0  for  gw  =  .5 
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Fig.  2f  Reynolds  Analogy  Factor  vs.  p  for  g  =  .5 
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Fig.  2g  Reynolds  Analogy  Factor  vs  P  for  g  =1.0 
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Fig.  2i  Reynolds  Analogy  Factor  vs.  (3  for  g  =1.0 
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Fig.  2j   Reynolds  Analogy  Factor  vs.  £  for  g  =1.1 
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Fig.  21  Reynolds  Analogy  Factor  vs.  £  for  g  =1.1 
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